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I. INTRODUCTION 

In recent years, there have been many theoretical stud- 
ies of geometrically frustrated systems^. These are sys- 
tems in which not all of the spin bonds can be satisfied si- 
multaneously due to the connectivity of the lattice. The 
frustration can lead to unconventional magnetic ordering, 
or to a complete absence of long-range order. 

Of particular interest are nearest-neighbor antifer- 
romagnets on bisimplej^ lattices composed of cor- 
ner sharing simplexes^. Examples include the two- 
dimensional kagome and the three-dimensional garnet 
lattices, each composed of triangular simplcxcs, and lat- 
tices composed of corner sharing "tctrahedra" : the two- 
dimensional checkerboard, the layered SCGOi, and the 
three-dimensional pyrochlore. 

The pyrochlore lattice, in which the centers of the 
tetrahedra form a diamond lattice^ (see Fig. 0), is of 
interest because it is realized in many experimental sys- 
tems, in both A2B2O7 oxides and in B sites of AB2O4 
spinels^, and because, by the analysis of Ref. 0, bisim- 
plex lattices composed of tetrahedra are less susceptible 
to ordering than lattices with triangle simplexes. 




FIG. 1: The pyrochlore lattice 



We consider the semiclassical limit of the pyrochlore 
nearest neighbor Heisenberg Hamiltonian 

u X./.a-s,. (1.1) 

ij 

where J^j = 1 for nearest neighbors (ij). and Si is a spin 
of length S^> 1 at site i. The Hamiltonian can be recast 
in the form 

H = ^|L a | 2 + const., (1.2) 

a 

where L a = ^ i£Q S;, is an operator that lives on the di- 
amond lattice sites a (we reserve Greek indices for tetra- 
hedra - diamond lattice sites - and roman indices for py- 
rochlore sites). The result is that, classically, all states 
satisfying 

L Q =0. (1.3) 

for all tctrahedra a, are degenerate classical ground 
states. Eq. (|1.3|) has three (L„,L^,L^,) constraints per 
tetrahedron, i.e. 3/2 constraints per site. Since there 
are two degrees of freedom (spherical angles 8,<f>) per 
site, there remain an extensive number of degrees of free- 
dom, and a degeneracy that is exponential in an exten- 
sive quantity^. This type of massive continuous degen- 
eracy has been shown in the kagome lattice antiferro- 
magnet to be lifted due to thermal^ or (anharmonic) 
quantum^ fluctuations, a phenomenon known as order 
by disorderiH^. However, in the classical case on the 
pyrochlore model, thermal fluctuations are not believed 
to facilitate orderingii^. 

In the S 3> 1 quantum model, one of us has 
demonstrated^ that the classical degeneracy is not fully 
lifted by the lowest order (in l/S) non-interacting spin- 
wave theory, assuming a collinear spin arrangement. In 
this paper, we recover this result using the more rigor- 
ous Holstein-Primakoff transformation, and provide a de- 
tailed study of various aspects in the linear spin-wave 
theory. 
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Recently, there have been various works designed to 
search for a ground state of the pyrochlore in the large-5 
limit. These include work on two-dimensional analogues 
of the pyrochlor o 15 ! 16 (see also Sec. lVIB|) . on ground state 
selection due to lattice distortion and spin-orbit coupling 
in Vanadium spinels^i, and on the large-A mean field 
theory^. Another body of work is on the closely related 
problem of the ordering in the pyrochlore in the presence 
of a magnetic field, which induces a collinear pattern^ 
(see also Sec. lVlA|) . 

The rest of this paper is organized as follows: In Sec. HP 
we derive the large-5 expansion of the Hamiltonian (|1.1|1 

H = TL C \ + Hharm + W qU art + • ■ ■ , (1.4) 

where H c \ is the classical Hamiltonian , of order 5 2 , 
Wharm ^ s the harmonic, non-interacting spin- wave, Hamil- 
tonian (of order 5) and H qua ,vt is the next order (0(1)) 
interaction termiS. We diagonalize the Harmonic part of 
the expansion, assuming fluctuation around a collinear 
classical ground state, where each state can be parame- 
terized by an Ising variable rji = ±1 on each lattice site, 
such that Si = r/iZ. 

In Sec. IIIII wc study the properties of the various spin- 
wave modes: zero modes that do not contribute to the 
zero-point energy, non-zero modes that can be expressed 
entirely in terms of the diamond lattice sites, and diver- 
gent zero modes that carry divergent fluctuations. We 
use this formalism to demonstrate a key result of Ref . Q: 
that collinear classical ground states related by Ising 
gauge-like transformations are exactly degenerate. 

In Sec. IIVI we consider the zero-point energy of the 
Harmonic fluctuations and derive an effective Hamilto- 
nian^, where we parameterize the energy only in terms 
of flux variables through each loop in the diamond lattice 
(where the pyrochlore spins sit bond centers). 

Kim = K+ A' 6 $ 6 + A 8 $ 8 + . . . , (1.5) 

where K n arc numerical coefficients that we evaluate us- 
ing a real-space loop expansion, and are sums of prod- 
ucts of the Ising variables r}% around loops of length n. 
We numerically evaluate the zero-point energy of a large 
number of collinear classical ground states and find that 
the effective Hamiltonian does a good job of evaluating 
the energy with just a few terms. We find a family of 
ground states and find an upper limit to the entropy, us- 
ing a correspondence between gauge-like transformations 
and divergent zero modes. 

In Scc.0wc demonstrate that the collinear states have 
lower energy than non-collinear states obtained by rotat- 
ing loops of spins out of collinearity, thereby showing that 
our assumption of collinearity was justified. In Sec. IVII 
we apply the loop expansion to some closely related mod- 
els including the case of non-zero magnetization plateaus 
and other lattices. We find that the effective Hamilto- 
nian approach is useful in predicting the ground states in 
many cases. In Sec. IVIll we present our conclusions. 



II. SEMICLASSICAL SPIN- WAVES 

In this section, we consider the effect of quantum fluc- 
tuations on the pyrochlore Heiscnberg model, in the semi- 
classical, S > 1 limit. In III Al wc perform a Holstein 
Primakoff transformation to expand the Hamiltonian in 
powers of 1/5. In 111 HI we focus on collinear classical 
ground states, Next, in Sec. Ill CI we diagonalize the 
collinear harmonic Hamiltonian to find the spin-wave 
modes. 



A. Large-,? expansion 

We start from a given ordered classical state, where 
the spin directions are parameterized by angles (fii,fa), 
such that the classical spin direction is 

n, = (cos 9i , sin 9i cos fa , sin Oi sin fa) . (2.1) 

We expand around this state, in powers of 1/5, to ac- 
count for quantum fluctuations, implicitly assuming here 
that the quantum fluctuations are small and do not de- 
stroy the local collinear order. Upon rotation to local 
axes (x,y,z), such that the classical spins are in the z 
direction (parallel to fij), we apply the usual Holstein- 
Primakoff transformation. Note that the choice of direc- 
tions x and y is arbitrary, and in the following we shall 
take y to be perpendicular to the z axis. One defines 
boson operators a,;, aj such that 

Sf = 5 — a\ai , 

5+ = S* + iS^ = ^2S-a\ai a t « 725a, , 

S- = S rc - ^5' 5 = a\yj2S-alai « V25aJ . (2.2) 

These operators satisfy the canonical bosonic commuta- 
tion relations 

[o», at] = Sij , [a»,Oj]=0, [aj,at] = 0. (2.3) 

We obtain the Hamiltonian l|1.4|l . where the leading term 
is the classical Hamiltonian 

n c i = 5 2 J « fi * • fi J » ( 2 - 4 ) 

y 

which is equivalent of Eq. Hl.lfl . and whose degenerate 
ground states satisfy 

J2 fi » = ■ ( 2 - 5 ) 

for all tetrahedra a. Due to the large classical degeneracy, 
we must go on to the leading order, harmonic, quantum 
correction, in order to search for a ground state, while as- 
suming that Eq. I|2.5|l is satisfied, i.e., that we are expand- 
ing around a classical ground state. The linear spin- wave 
energy E'harm was calculated by one of usi^, in previous 
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work, using classical equations of motion. The results of 
that work were that the classical degeneracy is not fully 
lifted by quantum fluctuations, to harmonic order, and 
that the remaining degeneracy is associated with a gauge- 
like symmetry. In this paper we justify these results using 
the more rigorous Holstein-Primakoff approach, which al- 
lows us to gain a better analytic understanding of the de- 
generacy, as well as perform numerical diagonalization. 
Furthermore, the Holstein-Primakoff transformation al- 
lows for a controlled expansion in powers of 1/5, includ- 
ing anharmonic order. In following work, we intend to 
show that the interacting spin-wave theory breaks the 
harmonic degeneracy^. 

We find it convenient to change variables to spin devi- 
ation operators 



n/_(a-aT), (2.6) 



^ = \/f(a + at) ; 



where at = (a j , a\, . . . , ajy ) is a vector of operators of 
length N s (the number of sites). We shall, from now on, 
reserve boldface notation for such vectors and matrices. 
These operators satisfy the commutation relations 



\a x ,<T v ] 



iSt 



\<J ,<T 



[cr v ,(7 y ] 



0. 



(2.7) 



The harmonic Hamiltonian can now be written in matrix 
notation 



((<7*)t,(<r^ 



R* P 



SN S , 

. (2 : 8) 

where the block matrixes, with respect to lattice site in- 
dex are 

Pij = -y- cos 6i sin 4>ij , 



Rfj = 5ij + (sin 6i sin 6j + cos 6{ cos 9j cos 4>ij ) 



Iff- 



(2.9) 



where we defined 4>ij = 



<j>j . These matrices depend on 
our arbitrary choice for the local transverse directions x 
and y, and can therefore not be expressed solely in terms 
of the spin direction n^. Note however, that in the case of 
coplanar spins, one can take, with no loss of generality, 
4>i = for all sites, i.e. spins in the (x, z) plane, and find 
that Pij = and that the matrix in Eq. I|2.8|l is block 
diagonal. 



B. Harmonic Hamiltonian for Collinear states 



is incompatible with collinearity. However, in the case 
of the pyrochlore lattice, or any other lattice composed 
of corner sharing tetrahedra, there is an abundance of 
collinear classical ground states. We will henceforth as- 
sume a collinear classical ground state, i.e. that each site 
is labeled by an Ising variable 7j,- = ±l, such that Si = r\{Z. 
We will return to the more general, non-collinear case in 
Sec^to justify this assumption a posteriori. 

In the collinear case (</>j = 0, 0i G {0, ir}), our definitions 
of the local axes for site i. are such that 



T]iX. 



V = V; 



(2.10) 



We may restore the x-y symmetry of the problem, by 
transforming the spin deviation operators back to the 
regular axes by the reflection operation 



r/rr 



, (2.11) 



where 



fm. o o 

?7 2 

773 

V 



2^ 

o 



(2.12) 



>}N 



The transformation (|2.11|l amounts to a reflection with 
respect to the y-z plane and makes the diagonal blocks 
in Eq. (|2.8|) equal to each other: 



H 
H 



(T-STrH. 



(2.13) 



where 



(2.14) 



H = 1 + — = -WtW, 

2 2 

J is the matrix whose elements are Jy and W Q i is an 
N s /2xN s matrix that takes the value 1 if i G a, and 
otherwise. 

The transformation of Eq. (|2.11() was chosen to ex- 
plicitly show the symmetry between the x and y axes in 
the collinear case. It also causes the Hamiltonian (|2.13|) 
to appear independent of {rji}. However, the particular 
collinear configuration does enter the calculation via the 
commutation relations 



r x rrV\ 



iSr] 



(2.15) 



Therefore, the equations of motion that govern the spin- 
waves do depend on the classical ground state configura- 
tion. 



The preceding derivation (Eqs. (|l.H2.<j[0 is valid for 
any lattice composed of corner sharing simplexes. One 
can argue on general grounds^*i2i22iS that the spin- wave 
energy has local minima for classically collinear states. 
In lattices that arc composed of corner-sharing trian- 
gles, such as the kagome, the constraint (Eqs. (|1.3I2.5|) '1 



C. Diagonalization of harmonic Hamiltonian 

Next we would like to Bogoliubov-diagonalize the 
Hamiltonian of Eq. I|2.13|l . so that we can study its eigen- 
modes and zero-point energy. As a motivation, we write 
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the equations of motion 
da 



dt 



-i[a,Hh a 



2S 



J7H 
r]H 



(2.16) 



These are the quantum equivalent of the classical equa- 
tions of motion derived in Ref^4. Upon Fourier trans- 
forming with respect to time, and squaring the matrix, 
we obtain 



\2S ) ° - { (r,H) 2 



(2.17) 



The spin-wave modes are therefore eigenvectors of the 
matrix (77H) 2 , with eigenvalues A 2 n = (w m /2S) 2 . If the 
non-hermitian matrix 77H is diagonalizable then its eigen- 
vectors {v m } are the spin-wave modes with eigenvalues 
{±|A m |}. Here v m is a vector of length N s , whose in- 
dices are site numbers. Note that although we refer to 
the eigenvectors of nil as the spin-wave modes, strictly 
speaking, the quantum mechanical modes are pairs of 
conjugate operators a x =w m , a y = isgn(A m )v m . 

Since {H 1 / 2 v m } are eigenvectors of the Hermitian ma- 
trix H 1 / 2 r/H 1 / 2 , the complete basis of eigenvectors {v m } 
can be "orthogonalized" by 

(v m ,T/v„) = v^rjvn = -i-v^Hvn = c m 8„^ n . (2.18) 

We shall henceforth refer to the operation (v m ,?7V n ) 
as the inner product of modes v m and v„. Note that 
Cm = (v m , fl^m) is not really a norm, since it can be zero 
or, if A m < 0, negative. The Bogoliubov diagonalization 
involves transforming to boson operators 



b m = 



1 



((jJVmjV 1 + isgn(c m )(j7v m ) t cr' y ) , 



V2S\Cn-, 

(2.19) 

with canonical boson commutation relations, to obtain 



Wharm — ^ ' 



m u m u ra 



bm + - 



SN S , 



with the zero-point energy 



Eharm — r W m — S I A^ 



SN S 



(2.20) 



(2.21) 



The fluctuations of cr x ' v are now easy to calculate from 
the boson modes 



S 



{a x (a x )^) = {a y (a^) = ^ VrnV t mj 



0. 



(2.22) 



Note that this is a matrix equation, where a x l y {a x l y )^ 
arc N s x N s matrices. 



III. SPIN- WAVE MODES 

We now examine the spin-wave modes that we found 
in Sec. Ill Cl to study their properties. In IIII Al we clas- 
sify them as zero modes, that do not contribute to the 
zero-point energy, and non-zero modes, that can all be 
expressed in terms of the diamond lattice formed from the 
centers of tetrahedra. Within the zero modes, we find, 

rV3% 



in IIII Bl that a number of modes (proportional to N s ) 
have divergent fluctuations. In Sec. IIII 01 we discuss the 
energy band structure and identify special singular lines 
in the Brillouin zone. 



In IIII Dl we show that the energy 
of any non-zero mode is invariant under an group of Z2 
gauge-like transformations. 



A. Zero modes and Non-zero modes 

Many of the eigenmodes of the Hamiltonian (|2.14|) , are 
zero modes, i.e. modes that are associated a frequency 
u> m = 0. These are the modes v 2 that satisfy 



(Wv z ) tt 



v z («) = for all tetrahedra a . (3.1) 



Since there are N s /2 tetrahedra, and thus N s /2 con- 
straints, we can expect as many as half of the spin-wave 
frequencies {w m } to be gencrically zero. In order to prove 
that indeed half of the eigenmodes of H are zero modes 
we first note that the subspace of zero modes is spanned 
by modes that alternate along loops, as in Fig. [2 i.e., for 
a loop £, up to a normalization factor 



(— l) ni i site number rij in loop £ , . 

iic ■ [6 - Z) 



In order for these loop modes to indeed be zero modes, 
two consecutive bonds in the loop cannot be from the 
same tetrahedron, i.e. the sites along each loop are cen- 
ters of bonds in a diamond lattice loop. From now on, we 
will use the term loop only for lines in the lattice satisfy- 
ing this constraint. The loop zero modes can in turn be 
written as a linear combination of hexagon modes only, of 
which there are N s (see Fig- E^) - Hexagons are the short- 
est loops in this lattice. Since there is a linear dependence 
between the four hexagons in a big super-tetrahedron (see 
Fig- El 3 ); then only N s /2 of the hexagon modes are lin- 
early independent 2 ^. Therefore a basis of zero modes of 
H would consist of half of the hexagon modes, and the 
matrix (rjH) 2 has at least N s /2 zero modes. Note that 
since these modes are zero modes of the matrix H, they 
are independent of the particular spin arrangement. We 
refer to these N s /2 modes as the generic zero modes, to 
distinguish them from other zero modes of (T7H) 2 that 
are not zero modes of H, that we shall discuss in the 
next section. All collinear ground states have the same 
generic zero modes. 

Zero modes mean that the quadratic correction to the 
classical energy due to small deviations from collinear or- 
der vanishes. It is interesting to note, that the zero modes 
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FIG. 2: (Color Online) (a) Example of a zero mode (in (001) pro- 
jection) alternating along a loop. Any such loop can be expressed as 
a linear combination of N s hexagon modes. Furthermore, there is a 
linear dependence between the hexagons in each super-tetrahedron, 
as depicted in (b). The additional constraints reduce the number 
of independent modes to N a /2. 

associated with loops that have alternating classical spins 
rji, are modes associated with the transformations under 
which 7i c i is invariant. In fact, deviations that rotate 
spins along these loops by any angle 6 take a collinear 
state to a non-collinear classical ground state, and when 
9 = 180°, the rotation takes one collinear classical ground 
state to another. These modes are completely analogous 
to the so-called weathervane modes in the kagomc lat- 
tice^. We will use the properties of these rotations in 
Sec. II V i^l to generate a large number of collinear classi- 
cal ground state, and in Sec. to study the dependence 
of the zero-point energy on deviations out of collinearity. 

Thus, the abundance of spin- wave zero modes is a re- 
flection of the macroscopic continuous classical degen- 
eracy. On the experimental side, the properties of the 
spinel material ZnC^Oi, at temperatures just higher 
than the phase transition into an ordered state, have 
been shown to be dominated by so-called local soft modes, 
which are spin- wave modes that have zero or infinitesimal 
frequency^. 

Whereas the zero modes are associated with the large 
classical degeneracy, they do not contribute to the quan- 
tum zero-point energy (|2.21|l . We now consider harmonic 
modes that have non-zero frequency. For any such mode 
v nz , the diamond lattice vector u = Wv„ 2 is an eigen- 
vector of (W?7Wt) 2 with the same eigenvalue 

(WryW^u = W(rjW f W) 2 v nz = A 2 u (3.3) 

In this fashion, we can get rid of N s /2 generic zero modes 
by projecting to a space that lives on the diamond lattice 
sites only, and is orthogonal to all generic zero modes. As 
already noted, since the generic zero modes are the same 
for all states, we can consistently disregard them and 
limit ourselves solely to diamond lattice eigenmodes of 
Eq. H3.3fl . Since the matrix Wr/W^ is symmetric, it is 
always diagonalizable and its eigenvectors are orthogonal 
in the standard sense. 

(3.4) 

We refer to these remaining N s /2 spin- wave modes, that 
can be viewed as diamond lattice modes, as the ordinary 
modes. Although the ordinary modes do not gcncrically 



have zero frequency, we may find that for a given classical 
ground state, some of them are zero modes. We will find, 
in the following section, that these are modes that have 
divergent fluctuations. 

B. Divergent modes 

As is apparent in Eq. I|2.22[l . divergent harmonic fluc- 
tuations occur whenever a certain mode satisfies 

Cd = (v d , f]v d ) = . (3.5) 

This can be shown to occur if and only if 77H is not 
diagonalizable, i.e., when the Jordan form of the matrix 
has a block of order 2. 

V Uv d = 0, 

r)Hw d = v rf . (3.6) 

We call the mode w<j, which is a zero mode of (r/H) 2 , 
but not an eigenmode of r/H, a divergent zero mode. As 
we noted before, a more accurate description would be to 
refer to the pair of conjugate operators <r x ~Vd, cr v ^Wd 
as the quantum mechanical divergent mode. From this, 
using Eq. (|2.14|) , we find that any divergent zero mode is 
related to a zero mode of Eq. I|3.3|l , which is defined on 
the diamond lattice 

u d = Ww d . (3.7) 

We call Ud a diamond lattice zero mode. Diamond lattice 
zero modes can be constructed so that they are restricted 
to only one of the two (FCC) sublattices. For example, 
one can construct an even divergent mode u^, i.e. 

u d (a) £ {0, ±1} for even a , Ud(ct) = for odd a , (3.8) 

such that for each of the adjacent odd sites /3 

u d(l)Vi(i,i3) = 0. (3.9) 

7: (7/3) 

Here 1(7,/?) is the pyrochlore site that belongs to both 
7 and (3. The mode is constructed in the following way: 
in addition to a, one other, freely chosen, (even) near- 
est neighbor 7 of /3 must be in the support of u^. If 
the bond connecting a and 7 is satisfied (antiferromag- 
netic) 1*^(7) = Ud{a). Otherwise, Ud{^) = —Ud{a) (see 
Fig. |2l • We repeat this until the mode span the sys- 
tem size, in the percolation sense, and all odd tetrahedra 
satisfy Eq. (|3~§|) . We call this an even divergent mode, 
and odd divergent modes arc constructed similarly from 
the odd sublattice. Even and odd divergent modes are 
linearly independent from each other, and any divergent 
modes can be separated to a sum of an even divergent 
mode and an odd divergent mode, and since odd and even 
modes are linearly independent from each other, one can 
construct a real-space basis for divergent modes from the 
singlc-sublattice divergent modes. 
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FIG. 3: (Color online) (a) Example of a portion of a (001) slice 
out of an even divergent mode, in real space (dashed lines connect 
to adjacent slices). Light (dark) colored bonds represent satisfied 
(unsatisfied) antiferromagnetic bonds, and an even tetrahedron a is 
marked "±" for u^(a) = ±1. Each participating tetrahedron must 
be have four others so that Eq. 1^.91 is satisfied for all (even and 
odd) tctrahcdra. Thus two even tetrahedra a, ■y have (a) = Ud (7) 
if they are connected by a satisfied bond, and u^ict) = —Udi'i) 
otherwise, (b) Example of a planar (xy) divergent mode, i.e. a 
mode that is bounded along the z axis (in this case, to one slice). 
Divergent modes can be bounded at most, along one of the major 
axes. 



If we examine the rules of constructing real space diver- 
gent modes, we find that their support (on the diamond 
lattice) can either be unbounded in space, or bounded 
along only one of the major axes, as in Fig. Eh- This 
can be easily proved: suppose that an even divergent 
mode is bounded in the x direction. There is an even 
diamond lattice site a that is in the support of this di- 
vergent mode, and its x value takes the maximum pos- 
sible value x a Note that in a slight abuse of notation, 
here (x,y,z) are spatial coordinates, while in Scc.lTTlthcv 
referred to spin directions. The site a has an odd neigh- 
bor at (x, y, z) = (x Q + l/4, y Q + l/4, z a — 1/4) (where the 
lattice constant of the underlying cubic lattice is taken 
to be 1), that has two (even) neighbors at x > x a (and 
thus cannot be in the support of the mode), and two 
neighbors at x~x a . One of these is a, and by the rules 
of the construction the other, at (x a , y a + l/2, z a ~ 1/2), 
must also be in the support of the mode. Continuing 
this reasoning, we will find that the site at coordinates 
{x a , y a + l, z a — 1) is also in the support of the divergent 
mode, and by induction the support of this mode is un- 
bounded in the +y and —z directions. Similarly, we can 
show that the support for this mode is unbounded in the 
— y axis and the +z directions. 



C. Magnon energy band structure 

Our discussion so far has focused on spin waves in real 
space, and has therefore been applicable to any (even 
non-periodic) classical collinear ground state. However, 
in order to perform numerical calculation (in Sec. II V Bj , 
one must consider periodic spin configurations, with a 
(possibly large) magnetic unit cell. We refer to the lat- 
tice composed of the centers of the magnetic unit cells 
as the magnetic lattice. We can Fourier transform the 



Hamiltonian l|2.13[) . using 

VWi^ q 

<?< = -^V^( r+i '), (3.10) 
q v^V 

where r is a magnetic lattice vector, / is a sublattice in- 
dex, corresponding to a basis vector A; , Nm is the num- 
ber of magnetic lattice sites, and q is a Brillouin zone 
vector. The elements of the transformed Hamiltonian 
matrix are 

JWq) = «im+^5> <q ' ftm ' ( 3 - n ) 

where the sum is over all nearest neighbor vectors & m 
connecting the sublattices / and m. Upon diagonaliza- 
tion, we obtain that the number of bands in the Bril- 
louin zone is equal to the number of sites in the mag- 
netic unit cell, i.e. the number of sublattices. The bands 
can be classified as follows: half of the energy bands be- 
long to generic zero modes, which have vanishing energy 
throughout the Brillouin zone. These modes are the zero 
modes of TC and are identical for all collinear classical 
ground states. The other half of the spin-wave modes are 
the ordinary modes, which can be viewed as diamond lat- 
tice eigenmodes of Eq. <|3.3|) . Of these bands, the optical 
ordinary modes possess non-zero frequency throughout 
the Brillouin zone, and the acoustic ordinary modes have 
non-zero energy in most of the Brillouin zone, but vanish 
along the major axes in reciprocal space (see, for exam- 
ple, the solid lines in Fig. llOf) . We find that the number of 
acoustic zero modes does depend on the particular clas- 
sical ground state, and the zero modes in these bands are 
the divergent modes, i.e. modes with divergent fluctua- 
tions (cTqClIq). 

Why are the divergent modes restricted to divergence 
lines in q space? Since the divergent modes are bounded, 
at most, in one direction, the Fourier space basis is con- 
structed by taking linear combinations of bounded planar 
modes, which are localized modes along the (major) axis 
normal to the plane. We label the (100) modes in such a 
basis by {w^}, where x is a (real space) coordinate along 
the normal axis, and the index m reflects that there may 
be several types of plane modes in each direction [for 
example, there could be four different types of square 
divergent modes of the form shown in Fig ^p] . Fourier 
space divergent modes can thus be constructed by linear 
combinations 

w (,.,o,o) = ^ w x e < to x j (3 . 12) 

X 

and similarly for y, z. These linear combinations can be 
taken at any q. value along the normal axis, and the 
number of divergent modes along each axis is equal to 
the number of values that the index m can take. The 
conclusion we can draw from this will be important later 
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on (in Sec. IIV C|) : the rank of the divergent mode space 

is of order N^ 3 - 

If we look at the acoustic energy bands to which the 
divergent modes belong, moving away from the diver- 
gence lines, in q space, we find (Fig. [TUjl that the en- 
ergy increases linearly with , the component of q per- 
pendicular to the divergence line. The dispersion of the 
acoustic modes can be easily found analytically by solv- 
ing Eq. (|3.3[) for small deviations away from a planar 
divergent mode (e.g. the one depicted in Fig. El 3 ) ; with 
q± restricted to be within the plane. 

The singular spin fluctuations along lines in q space 
would produce sharp features in the structure factor 
5(q), that could be measured in elastic neutron diffrac- 
tion experiments. This is because the structure factor 
is proportional to the spin-spin correlation (Sq • S_ q ), 
where S -1 is the component of the spin transverse to the 
scattering wave vector q. To lowest order in S, only a x 
and a y contribute to the correlations, and one obtains 

S(q)oc£<a q <rI"q>. (3.13) 

(lm) 

Thus, the structure factor should have sharp features for 
any q along the major lattice axes. 

Furthermore, since the zero-point energy of the present 
harmonic theory will be shown to have degenerate ground 
states (see Sec. UTTPl and Ref. Q, anharmonic correc- 
tions to the harmonic energy determine the ground state 
selection^. It turns out that the divergent modes become 
decisive in calculating the anharmonic energy 7i qU art- 
The anharmonic spin-wave interaction would also serve 
to cut off the singularity of the fluctuations. 

In Sec. IIV CI we shall show that the divergent modes 
also provide a useful basis for constructing and counting 
gauge-like transformations that relate the various degen- 
erate ground states. 

D. Gauge-like symmetry 

Upon examination of Eq. (|3.3|l . it becomes apparent 
that the harmonic energy -Eharm, of Eq. I|2.21[) . is invari- 
ant under a Zi gauge-like transformation that changes 
the sign of some tetrahedra spin deviations^ 

m T(a) T ((3) m , (3.14) 

where r(a),r(/3) € ±1, and a, f3 are the two tetrahedra 
that share site i. While this is an exact gauge symmetry 
of the projected (diamond lattice) Hamiltonian (related 
to 13. 31 ). it is not a physical gauge invariance, since the 
transformation must be carried out in a way that con- 
serves the tetrahedron rule, i.e. does not take us out of 
the classical ground state manifold. Furthermore, these 
transformations relate physically distinct states. We will 
henceforth refer by transformation to "allowed" gauge- 
like transformations that conserve the classical tetrahe- 
dron rule. 



Fig. 0] shows an example of a transformation that in- 
cludes flipping some of the even tetrahedra only (even 
transformation). There are two transformations that can 
each be viewed either as an odd or as an even transforma- 
tion: the identity transformation (r(a) = 1 for any a) the 
global spin flip (r(a) = — 1 for any even a and t(q) = +1 
for any odd a, or vice- versa). 

For any even transformation r cvon , we can define its 
reverse (even) transformation by 

Teven(a) = -T cvcn (a) , V even a , (3.15) 

and similarly for odd transformations. Examining 
Eq. H3.14(l we find that there are actually two ways of 
expressing any transformation as an product of an even 
transformation and an odd transformation because 

^"even ® ^"odd = 7"even & 7"odd j 

(3.16) 

Even and odd transformations commute, in the sense 
that applying an even (odd) transformation does not af- 
fect the set of allowed odd (even) transformations. Thus 
we find that for any reference state, the number of pos- 
sible gauge transformations A/"<l re ^ satisfies 

A/£ cf) = jA® x M^ d f , (3.17) 

where A4ven and Af^P are the number of even and odd 
transformations for the reference state, respectively, in- 
cluding both the identity transformation and the overall 
spin flip. 




FIG. 4: (Color online) (001) projection of a portion of an even 
transformation. The shaded tetrahedra, on the even sublattice of 
the diamond lattice, are being flipped. Each flipped tetrahedron 
must be connected by 4 satisfied (highlighted) bonds to neighboring 
tetrahedra in the same sublattice. Up (down) spins are signified by 
closed (open) circles. Dashed lines connect one layer to another in 
the three-dimensional structure. 

If we take a particular classical ground state, and ap- 
ply a gauge transformation (as depicted in Fig. 0} to it, 
it is easy to see graphically what would happen to the 
divergent modes (as shown in Fig. |3J): If a tetrahedron 
marked by "+" in Fig[3]overlaps the support of the gauge 
transformation, it turns into a "— " and vise versa. More 
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formally, if we start from a given state, in which there 
is a divergent mode u^, a transformation r (such that 
r a =±1) results in the new state with a divergent mode 

u(a) — > r a u(a) , (3.18) 

for each a. Otherwise, the number and spatial support 
of the divergent modes is gauge- invariant. On the other 
hand, one finds that states that are not related by gauge 
transformations, and therefore have different energies, 
generally have a different number of divergent modes. 



IV. ZERO-POINT ENERGY AND EFFECTIVE 
HAMILTONIAN 

In this section, we study the zero-point energy of the 
harmonic Hamiltonian (|2.13|l . First, in Sec. IIV Al we 
write the energy as an expansion in real-space paths, 
and use this expansion to produce an effective Hamil- 
tonian in terms of Ising like fluxes through the diamond 
lattice loops. Next, in Sec. IIV Bl we numerically diago- 
nalize the Hamiltonian to calculate the zero-point energy 
for various classical ground states. We compare the re- 
sults to the effective Hamiltonian and find that they agree 
well, and predict the same ground state family. Finally, 
in Sec. IIV CI we find an upper bound for the number of 
harmonic ground states using a correspondence between 
gauge transformations and divergent modes. 



A. Effective Hamiltonian 

As the zero-point energy is gauge invariant, it can 
only depend on the gauge invariant combinations of {rji\, 
which are products of the Ising variables around loops. 
These can be viewed as Z 2 flux through the plaquettes of 
the diamond lattice^. In this section, we find an effective 
Hamiltonian in terms of these new degrees of freedom, us- 
ing a real space analytic loop expansion in the spirit of 
the l oop expan sions used in the large- N calculations of 
fiefs. Il8l27l2sl 

The harmonic energy of Eqs. (|2~2T|) . (j2~T31 and (|2~Ti|) 
can be written as 

^harm^STr^/X 2 ) 1 / 2 -^, (4.1) 

where [i = Wt7W t is the N s /2 x N s /2 matrix, whose 
indices are diamond lattice sites 

fJ-ap = w txiWf3ir]i • (4.2) 

i 

For any collinear classical ground state, the diagonal 
part of /it vanishes, elements connecting diamond near- 
est neighbors are equal to ±1 all other elements are 0. 
Therefore, the non-diagonal elements of /x 2 connects be- 
tween the same diamond sublattice, i.e. between FCC 



nearest neighbors. 

( 4 a = /3 

{^ 2 )aj3 = { Vctp a :P next nearest neighbors , (4.3) 

I otherwise 

where rj a p = rjirjj and (ij) is the (pyrochlore) bond con- 
necting a and /3. Thus, we could formally Taylor-expand 
the square root in Eq. (|4.1|) about unity. In order to as- 
sure convergence of the expansion, as will be discussed 
later, we generalize this and expand about At, where A 
is an arbitrary disposable parameter. 

£ harm = STi{At + {^- At)) 1 ' 2 -SN S (4.4) 

,.2 

= 5v / A^C„Tr(^--41)" - SN S 
= SVA^C^t^hW^-SNs, 
with the coefficients 

Or, _ 3V' 

C = l, C n = (-l) n+l [ \ n j " ,forn>0. (4.5) 

on! 

Tr fi 2k is a sum over all of the diagonal terms of fi 2k , i.e. 
a sum over products of /x Q)3 along all of the closed paths 
of length 2k in the diamond lattice. Here a "closed path" 
is any walk on the diamond lattice that starts and ends 
at the same site. This expansion involves constant terms 
that are independent of the sign of ii a p as well as terms 
that do depend on particular spin configurations. For 
example, any path of length 2k, such that each step in 
one direction is later retraced backwards (a self- retracing 
path), will contribute 1 to Tr/i 2fe . On the other hand, 
paths involving loops on the lattice could contribute ei- 
ther + 1 or —1 depending on the spin directions. Thus, we 
can re-sum Eq l|4.4|l to obtain an effective Hamiltonian 

K% m = Eo+K^+K^+Y, #io, s $io, s +- ■ ■ , (4.6) 

s (10) 

where Kn, i^2Z,s are constants, which we calculate below, 
and <f>2i (&2i,s) arc sums over all loops of length 21 (and 
type s). The index s is to differentiate between different 
types of loops of length 21 that are not related to each 
other by lattice symmetries. In our case, we do not need 
the index s for loops of length 6 or 8, since there is just 
one type of each. On the other hand, there are three 
different type of loops of length 10, and therefore there 
arc three different 21 = 10 terms. Since we will explicitly 
deal with just the first three terms in Eq. Q4.6[l . we shall 
omit the index s from now on. In fact, in the approxi- 
mation that we present below, we will assume that all of 
the loops of any length 21 have the same coefficient K21. 

By our definition of loops, <I>2Z can be expressed ei- 
ther in terms of 21 diamond lattice sites along the loop 
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(ax, Q!2j • ■ ■ ! Of2fc)i or i n terms of 21 pyrochlore lattice sites 
(ii,i2, ■ ■ ■ , «2l) 

(ai,Q2,---,a2i) 

XI ViiVh • • • Vi* ■ ( 4 - 7 ) 

(ll ,l2,...,i2j) 

Note that in general, paths that we consider in the cal- 
culation (and that are not simple loops) should only be 
viewed as paths in the diamond lattice. 



1. Bethe lattice harmonic energy 

Before we evaluate the coefficients in Eq. (|4.(j[) for the 
diamond lattice, we shall, as an exercise, consider the 
simpler case of a coordination z = 4 Bethe lattice. In 
order for this problem to be analogous to ours, we assume 
that the number of sites is Nb=N s /2 and that each bond 
(a/3) in the lattice is assigned an Ising variable ^ Q( g = ±1. 
In this case, since each bond included in any closed path 
along the lattice is revisited an even number of times, 
and since //^ = 1 for any bond in the lattice, then each 
closed path of length 2k contributes 1 to Tr/x 2fe . Thus 
all bond configurations in the Bethe lattice would have 
the same energy. 

Calculating the Bethe lattice energy for a given path 
length 2k turns out to be a matter of enumerating the 
closed paths on the Bethe lattice, which can be done 
exactly using simple combinatorics (see App. [SJ. The 
sum (|4.4|) does not converge, as we consider longer and 
longer paths, for the trivial choice of A= 1, but converges 
well for A > 1.4 (see Fig.[SJ). In the thermodynamic limit, 
we obtain (see Eq. (|A4|0 



^harm (Bethe) = E = -0.5640SW S 



(4.8) 



This value was obtained from Eq. I|4.4|l . cutting it off at 
n = 30 and extrapolating to n — > oo (see Fig.EJ). 



2. Bethe lattice approximation for the constant term 
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FIG. 5: (Color online) The analytically calculated constant term 
Eo, in the Bethe lattice approximation as a function of the maxi- 
mum path length considered, for various values of A. We find that 
the sum converges for A>1.4. 
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FIG. 6: (Color online) Extrapolation of the energy calculated for 
paths of length n < 30 to n — » oo, using the linear dependence of 
-Eo on 1/n. The calculated energy is shown in open circles and the 
extrapolated inverse linear dependence in lines. We see that the 
extrapolated results is nearly independent of our choice of A> 1.4. 



What can we learn from the Bethe lattice exercise 
about the diamond lattice effective energy? It turns out 
that actually the coordination-4 Bethe lattice calcula- 
tion provides a very good approximation for the constant 
term Eq. There is a one-to-one correspondence between 
the Bethe lattice paths, and self-retracing diamond lat- 
tice paths, which arc the paths that contribute to Eq. 
Conversely, the product of \i a ^ along a path that goes 
around a loop depends on the particular classical ground 
state, and therefore contributes to the term in Eq. (|4.6|) 
corresponding to that loop and not to Eq. 

There is only one type of path that contributes to the 
constant term, but was omitted in the Bettc lattice ap- 
proximation: If a loop is repeated twice (or any even 



number of times), in the same direction, as in Fig. [7|d, 
it does contribute a constant term since each fi^p = 1 
for each link along the loop. A simple argument can be 
given to show that the constant terms involving repeated 
loops are negligible compared to the Bethe lattice terms. 
First, we realize that the contribution of paths involving 
repeated loops is exactly the same as the contribution of 
paths that have already been counted by the Bethe lattice 
enumeration: self-retracing paths that go around a loop, 
and then return, in the opposite direction (see Fig. UK)- 
We call these self-intersecting self-retracing paths since 
they include retraced diamond lattice loops. Since they 
still correspond to a Bethe lattice paths, they can be 
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FIG. 7: Repeated loops that contribute constant terms: (a) Cir- 
cled back and forth, and accounted for by the Bethe lattice enu- 
meration (self-intersecting paths), and (b) Repeated in the same 
direction, not accounted for in our calculation. 

readily enumerated. 

In order to get an idea for the number of such self- 
intersecting paths, consider the smallest loop, a hexagon. 
The number of paths of length 2k = 12 starting from 
a particular point is (by App. [XJ, / 6 = 195352, while 
the number of hexagons touching that point is only 12, 
accounting for 24 intersecting paths of length 2k < 12! 
If we now consider the paths of length 14 starting from 
the same point, then there are less than 100 involving 
a repeated hexagon, while there are nearly 2 million in 
total. Essentially, the number of self-intersecting paths, 
involving loops of length 21 is smaller by a factor greater 
than (z — I) 1 than the total number of paths of the same 
length. 

The argument is reinforced a posteriori, by enumer- 
ation of paths involving loops, that we do below. The 
lowest order correction to the constant term Eq (from 
Eq. H4.8fl due to repeated loops is of the same order as 
the coefficient for a loop of particular loop of length 12 
(circling twice around a hexagon) , which we find to be of 
order 10 -4 . 



3. One loop terms 

We now move on to calculate the coefficients K21 of the 
non-constant terms that involve simple loops. The pre- 
diction of Ref.^Jis that these terms decay with increased 
length of loops and therefore an effective Hamiltonian of 
the form i|4.6|) can be derived. 

Consider a particular loop of length 21. We try to enu- 
merate all of the closed paths of length 2k that involve 
this loop and no other loops, i.e. all of the terms propor- 
tional to Yli=i A^ctiQi+ii where a, are the sites along the 
loop (a 2 i+i=a 1 ). 

Since we allow for no additional loops, we assume that 
the path is a decorated loop, i.e. a loop with with self- 
retracing (Bethe lattice-like) paths emanating from some 
or all of the sites on it. In order for this description to be 
unique, we allow the self-retracing path emanating from 
site on along the loop to include site Q!i_i, but not site 
Oi+i. Thus, the first appearance of the bond (i, i + 1) is 
attributed to the loop, and any subsequent appearance 



must occur after going back from a site j > i, and is 
attributed to the self-retracing path belonging to j. 

App. [5] describes the practical aspects of this calcu- 
lation. See Fig. 1151 for a diagrammatic description of 
the paths we enumerate. The approximation neglects, as 
before, the contribution of repeated loops, which is neg- 
ligible, by the argument of Sec. II V A~2l This means that, 
within out approximation, all loops of length 21 have the 
same coefficient in Eq. (|4.6ll . Calculating the sums in 
Eq. (|4.4|) for n < 30 and extrapolating to n — » 00, we 
obtain the values A" 6 = 0.01365, K s = -0.00335". 

Thus we have evaluated the first three coefficients in 
the effective Hamiltonian <|4.6[) . and in the next section, 
we shall compare the effective energy to numerical diag- 
onalization. 



B. Numerical diagonalization 

In order to be able to test our predictions numerically, 
we first constructed a large number of classical ground 
states, using a path flipping algorithm — on a cubic unit 
cell of 128 sites, with periodic boundary conditions. In 
each step of this algorithm we randomly select a loop of 
alternating spins to obtain a new classical ground state. 
Considering the large classical degeneracy, we can con- 
struct a very large number of collincar classical ground 
states in this manner. In order to explore diverse re- 
gions of the Hilbert space, we started the algorithm with 
various different states that we constructed by hand. 

We have Fourier transformed the Hamiltonian l|2.13|l . 
with a magnetic unit cell of 128 sites, diagonalized for 
each q value, and calculated the harmonic zero-point en- 
ergy H2.21f) for fluctuation around a wide range of classi- 
cal ground states. We show the calculated -Eharm for 50 
sample states in Fig. Our calculations verified that 
indeed gauge-equivalent states always have the same en- 
ergy. 

We also show in Fig, |S] the effective energy E^ rm for 
the same 50 states, using Kq and Kg, calculated above. 
The effective Hamiltonian, even with just 3 terms, proves 
to do a remarkably good job of approximating the energy. 
The root-mean-squared (RMS) error for the 50 states 
shown in Fig. [SJ is 1 X 10~ 3 , and it can be attributed to 
higher order terms (21 > 8) in the expansion 1)4. 6fl . An 
independent numerical fit of the three constants Eq, Kq 
and Kg can lead to a slightly smaller RMS error, of the 
same order of magnitude, but the coefficients are strongly 
dependent on the set of states used in the fit. 

Based on the functional form of Eq. (|4.6|) . it has been 
speculated 1 ^, that the ground state manifold consists of 
all of the (gauge-equivalent) so-called ir-flux states (using 
the terminology of Ref. fl5|) , in which 

n ^ = -i, (4.9) 

for all hexagons. We indeed find that there is a family of 
exactly degenerate ground states, satisfying Eq. I|4.9|l . In 
Fig. we show some of the ground states. The smallest 
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FIG. 8: (Color online) Calculated zero-point energy for 50 ran- 
domly generated classical ground states, compared to the Ising ef- 
fective energy of Eq. 14.61 . with 21 < 8. State number 1 in the plot 
is the 7r-fiux state and state 50 is the zero-flux state, in which all of 
the terms <J>2! take the maximum possible value. Notice that the 
effective Hamiltonian is not as good at calculating the energies of 
these two extreme states, as it does of calculating the energies of 
other states. This is because neglected higher order terms in these 
states tend to add up in these states rather than cancel out. 



magnetic unit cell that we obtain for a ground state has 
16 spins, although if we consider bond variables rather 
than spins (as in Fig. |5J), the unit cell can be reduced 
to 8 sites. The highest energy, among collinear states, is 
obtained for the Q-flux states, for which the spin direc- 
tions have a positive product around each hexagon. Note 
that since the spin product around any loop can be writ- 
ten as a product of hexagon fluxes, all states satisfying 
Eq. I|4.9|) have the same loop terms <&2z, for any I, and 
are thus gauge-equivalent. In general, if two states have 
the same hexagon product for every hexagon, they are 
necessarily gauge-equivalent. 

To search for ground states, we developed a computer 
algorithm that randomly generates even or odd gauge 
transformations starting from a particular state, and a 
unit cell with periodic boundary conditions. We start 
by flipping a random even tetrahedron. In each subse- 
quent step, we find an odd tetrahedron that violates the 
tetrahedron rule (i.e. has non-zero sum), and randomly 
flip one of its (even) neighbors that can fix that viola- 
tion. This is repeated until there are no more violated 
tetrahedra. A similar algorithm is employed to find odd 
transformations. 

We performed an exhaustive search for ground states, 
satisfying Eq. (|4.9|) . in the 128 spin cubic magnetic unit 
cell of linear dimension L = 2. We started with a par- 
ticular 7r-flux ground state and randomly generated 10 
even and 10 7 odd gauge transformations and found 142 
unique transformations on each sublattice, resulting, by 
Eq. (joTTI) . in a total of 10082 distinct states. Only 
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FIG. 9: (Color online) (001) slices of some of the ground states 
with smallest magnetic unit cells. For clarity, we only show the 
bond types, which are (unlike the spin directions), identical for all 
(001) slices, for these particular states,. Ferromagnetic bonds are 
presented in dark and antiferromagnetic bonds are light. Dashed 
lines connect this slice to adjacent slices. 



24 of these are unique with respect to lattice and spin- 
flipping symmetries. Note that although the number of 
even gauge transformations and the number of odd gauge 
transformations turn out to be the same for the harmonic 
ground states, they need not be the same for other states. 
Fig|31 shows the states with smallest unit cells. By con- 
struction, all of these states are exactly degenerate to 
harmonic order in spin-wave theory. In future work, we 
shall explore the anharmonic selection among the har- 
monic ground states^i. 



C. Ground state entropy 

In the preceding section we have found that there is 
a large family of exactly degenerate ground states. Any 
two of these states are related by a gauge transformation 
of the type discussed in Sec. lIIIDl and therefore, in order 
to enumerate these states, we must find how many gauge 
transformations one can perform on a give n ground state, 
in an arbitrarily large system. In Ref. Il4l the number of 
ground states was speculated to be of order e L , where 
L is the linear dimension of the system size (as opposed 
to a classically extensive entropy). However, that was 
only rigorously shown to be a lower bound, by explicitly 
constructing a set of layered ground states in which each 
layer can independently be flipped. Here, we aim to find 
an upper bound for the number of ground states. 

Comparing Figs. 0] and [3J we see that there is a close 
relation between even (odd) transformations and even 
(odd) divergent modes. Any (diamond lattice) divergent 
mode that is constructed with satisfied antiferromagnetic 
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bonds only, describes a valid gauge transformation. This 
is reminiscent of the relation that we saw in Sec. 1111 Al 
between generic zero modes of 7iharm and transforma- 
tions that keep 7i c \ invariant. The relation between di- 
vergent modes and transformations, as illustrated above, 
and our knowledge of the divergent modes allow us to 
demonstrate that the entropy has an upper bound of or- 
der Lin L. 

Consider a particular reference harmonic ground state. 
We assume that all of the ground states are related to 
each other by gauge transformations. Thus, any ground 
state is related to the reference state by a transformation, 
that can be almost uniquely expressed as product of an 
even and an odd transformation. The number of ground 
states TVo is equal to the number of transformations pos- 
sible for the reference state Nq C1 ^ (see Eq. (13.171) ). 

M)=A/£ cf) . (4.10) 

We have seen in Sec. IIII CI that the number of indepen- 
dent divergent modes is of order L. In the 7r-flux har- 
monic ground states, we find that there are 24L indepen- 
dent divergent modes, where L is measured in units of 
the underlying cubic lattice, half of which are even and 
half odd. Consider an orthogonal basis of 12L real space 
planar even divergent modes, i.e. {u™ cn } are a set of 
unnormalized orthogonal vectors living on the even di- 
amond sublattice, satisfying Eqs. (|3.4|) and (|3.8() . Any 
even transformation is associated with a divergent mode 
u even a l so living on the even sublattice. We can write the 
transformation in terms of the basis 



even / _j 
m,—l 
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(4.11) 



Each of the planar divergent modes has a support of 
2L 2 diamond lattice sites, and since u™ en (a) is or ±1 
for any a (Eq. G£EJ>), we find that |u™ on | 2 = 2L 2 Since 
u ovon( a ) is or 1 for each a, the inner product u^, en -u™ on 
is also an integer satisfying 



(4.12) 



Therefore, each coefficient in Eq. I|4.11|l can take one of, 
at most, 4L 2 values, and since there are 12L basis vectors 
u™ en , the number of possible vectors u^, cn is, at most, 
{2L) 2AL . This is an upper bound on the number of even 
transformation, and similarly, of odd transformations, as 
well. By Eq. (|5T7|) . we find 



■M, < ^(2L) 



48L 



(4.13) 



The entropy is defined as lnA/" , and is, at most, of order 
L\nL. From Ref. [3 we know that Afo > 2 4 L, and the 
entropy is at least of order L. 

This same bound on the order of the multiplicity ap- 
plies to any family of gauge-equivalent Ising configura- 
tions, since we enumerated the possible gauge transfor- 
mations on any given reference state (not necessarily a 



7r-flux state), which implies that the upper bound to the 
number of states in any gauge-equivalent family is of the 
same order. However, while the order of magnitude of the 
multiplicities of all energy levels are the same, the coeffi- 
cients in front of LlnL differ, because the the number of 
independent divergent modes (which is always of order 
L) is generally not the same for different gauge families. 



V. NON-COLLINEAR SPINS 

We now move on to discuss the case of non-collinear 
classical ground states, aiming to show that the energy of 
any collinear ground state increases upon rotating some 
of the spins out of collinearity. 



A. Collinear states are extrema of Sharm 

Consider first the case of a coplanar rotation, i.e., with- 
out loss of generality, some of the spins are rotated such 
that the angle Bi is neither nor 7r, while <\>i remains 0. 
The elementary way of performing such a rotation is to 
rotate the spins in an alternating (in r\i) loop, by +9 and 
—9 in an alternating fashion, where 9 is a constant. 



9i = 

9r =7T 



(5.1) 



Carrying through the derivation of Eq. I|2.17|) . for the 
coplanar case, we find that the dynamical matrix ele- 
ments change 



//,///. //,/. ///. , -> cos(6».; - 9 k )H ik H kj 



(5.2) 



This means that in the expansion (|4.4(l , fi 2 changes from 
Eq. O to 

'4 a = /3 

{v?)u(} = \ cos 9 a p a , (3 next nearest neighbors . 
otherwise 

(5.3) 

Here 9 a p = 9i — 9j, where (ij) is the (unique) pyrochlorc 
bond a site i £ a and a site j G (3. 

To see how the zero-point energy changes with copla- 
nar deviations away from a collinear state, we take the 
derivative of the Hamiltonian l|2.8|) with respect to 9i 



dH ha 



39, 



= - Jy sin(6>i - 6> 3 >f of. 



(5.4) 



This is clearly a zero operator for any collinear state. 
Therefore (i) all collinear states {9i—9j £ {0, it}) are local 
extrema of £'harm({^})- (h) For small deviations of 6^ 
from collinearity, the deviation of the energy from the 
value is quadratic: 

£Wm({0i}) - ^harm(O) = C({(^ - 9j) 2 }) . (5.5) 
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expected to come from the divergent modes. We show 
below that the divergent modes' frequency is non-zero 
for coplanar states, and therefore the zero point energy 
increases upon rotation from a collinear to a coplanar 
state. 



B. Spin- wave modes upon deviation from 
collinearity 



0. 



-0.5 





(7u/8,qy,0) 



0.5 



FIG. 10: (Color online) Dispersion along a line in q space, for 
the collinear state shown in Fig. (solid lines), as well as the 
# = 7r/6 coplanar state (dashed lines), showing that the divergent 
zero mode gains non-zero energy along the (x,0,0) line, when the 
state is taken away from collinearity, while all other modes are 
virtually unchanged. Insert: Zero-point energy of the same state 
with the spins along the ferromagnetic x-y diagonal rotated by 
angle ±0. The minimum energy is for the collinear states. 



We have shown that the collinear states are local ex- 
trema of -Eharm, in the space of coplanar classical ground 
states. However, are they minima or maxima? In order 
to find out, we take the second derivative of the Hamil- 
tonian, evaluated at a collinear state 



d 2 H h; 



dOidOi 



coll 



JtfCOB^i-^Xof) 2 
k 

J ij (of) 2 - Sij Y J ik^i^k > ( 5 - 6 ) 
k 

— 6j) — > rjiTjj at the collinear state 
and used Eq. (|2.11|> to transform to the a x operators. 

Now we can evaluate the elements of the Hessian ma- 
trix M by taking the expectation value of this operator 



where we took cos( 



Mi. 



d 2 E ha 



dOiddj 



i 7 



'I = J 

otherwise 



• (5-7) 



Here we used Eq. <|2.22|1 and the properties of the eigen- 
modes to simplify the expressions. In order to show that 
the collinear states are local minima with respect to {0i\, 
we would have to show that M is positive definite. Al- 
though we have not been able to shown this analytically, 
we have found this to be true for a large number of peri- 
odic collinear states that we have constructed (where we 
introduce a cutoff to the singularity of the fluctuations). 

In order to examine the angle dependence in a gen- 
eral case, we note that the matrix elements My are 
dominated by the divergent fluctuations of the divergent 
modes. Therefore, the leading order change in -Eharm is 



In order to better understand the origin of the 
quadratic (in 9) energy change, we examine the eigen- 
modes of the new dynamical matrix (|5.2|) . All of the 
N s /2 generic zero modes of the collinear dynamical ma- 
trix are zero modes of R y = H in the coplanar case 
(using the local frame chosen in Sec. Ill Al see Eq. i|2.8|0 
and therefore remain zero modes for any coplanar spin 
arrangement. One finds, however, that these generic zero 
modes acquire divergent fluctuations in coplanar states, 
because of the difference in stiffness between in-plane and 
out-of-plane fluctuations. This is similar to the case of 
the kagome lattice, where all of the zero modes of copla- 
nar classical ground states have divergent fluctuation. 

On the other hand, the divergent zero modes of the 
collinear states, become (non-divergent) nonzero modes 
when the spins are rotated out of collinearity. If a loop in 
a collinear state is rotated, as in (|5.1|) . by ±0, we find that 
the divergent modes' frequency increases linearly with 9. 
This rise in the divergent modes' zero-point energy is 
the reason that each collinear classical ground state has 
lower energy than nearby coplanar states. However, af- 
ter integration over the Brillouin zone, the rise in total 
zero-point energy is quadratic in 9. The inset in Fig llUI 
shows a numerical example of this for the state in Fig.|51i, 
where all of the spins along the ferromagnetic line, in the 
(110) direction, are rotated by a small angle away from 
the axis of the other spins. Clearly, the lowest energy con- 
figurations arc collinear. In Fig.^| we see that the most 
significant difference in the dispersion between collinear 
and non-collinear states, is the gap formed in some of 
the zero modes along the divergence lines. The change 
in energy of the optical modes is negligible compare to 
this gap, for all cases that we looked at. 

Looking at further rotation of spins out of a coplanar 
arrangement, i.e. rotation of a loop or line by angles 
±0, one finds that some of the zero modes gain non- 
zero frequency proportional to \<f>\, as for the kagome 
lattic e) 30 ! 31 , and the energy increase is 0({<fii}). 

While we have shown that collinear states are local 
minima of the energy landscape, we have not ruled out 
the possibility that a non-collinear state would be a local 
(or even global) minimum. Although we do not believe 
this to be the case, further work would be required to 
prove so. 
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VI. 



EFFECTIVE HAMILTONIAN FOR 
RELATED MODELS 



The loop expansion of Sec. IIV Al can be easily be 
adapted to study similar models that support collincar 
ground state including the pyrochlorc Heisenberg model 
with a large applied magnetic field and the Heisenberg 
model on closely related lattices, as we shall show in 
Sees. I VI Al and IVIBI respectively. 

For the purpose of determining the ground state mani- 
fold, we find it convenient to recast the effective Hamilto- 
nian of Eq. (|4.6|) as an Ising model in the complementary 
lattice. This is a lattice composed of the centers of the 
shortest loops, or plaquettes, in the original lattice. In 
the pyrochlorc lattice, the complementary lattice sites 
are the centers of the hexagons, and form another py- 
rochlore lattice. To each complementary lattice site a, 
we assign an Ising spin f] a , equal to the product of the 
direct lattice sites around the corresponding plaquette. 



Urn, 



(6.1) 



where the product is on all sites in the (direct lattice) 
plaquette a. Since any loop product can be written as 
a product of spin products around plaquettes, the terms 
$2Z is Eq. (|4.6|1 are now represented by simpler (at least 
for small I) expressions, in terms of the complementary 
lattice spins 

C = £ + SB ]T Va + SJ VaVb 

a (ab) 

+sj' W + S ^ 

VaVbVc + ■ ■ (6,-2) 

((ab)) Aabc 

where (■ • • ) and ((• • ■)) represent nearest neighbors and 
next-nearest neighbors on the complementary lattice, re- 
spectively, and the A sum is over 3-spin plaquettes. Com- 
paring to Eq. I|4.6|) . we can identify, for the pyrochlorc, 
£ Q = Eq, B = K 6 /S, J = K S /2S, J' = K w /S, J 3 = 0. 
The 1/2 factor in J stems from the fact that each 8- 
loop in the pyrochlorc has two different representations 
as a product of two hexagons, and J-$ vanishes because 
in pyrochlorc complementary lattice the product of three 
spins of a tetrahedron is equal to the fourth spin in the 
same tetrahedron. This is due to the dependence, in the 
direct lattice, between the four hexagons in one super- 
tetrahedron (see Fig. . Thus the J3 term is already 
accounted for in the "field" B term. Note that, in the py- 
rochlore, there are three different types of loops of length 
10, and thus there will be 3 different terms in Eq. (|6.2|1 
that have a coefficient equal to J' , within the decorated 
loop approximation. 

Writing the effective Hamiltonian in terms of the com- 
plementary lattice spins is manifestly gauge invariant, 
since the complementary spins f] a are not modified by 
gauge transformations. In the pyrochlore Heisenberg 
model, we found in Sec. II VI a ferromagnetic nearest neigh- 
bor interaction J , with a positive field B, resulting in 



a uniform complementary lattice ground state in which 
fj a = — 1 for all a (the 7r-nux state) . This unique comple- 
mentary lattice state corresponds to the family of direct 
lattice ground states satisfying Eq. (|4.'jp . 



A. Non-zero magnetic field 

We now consider what happens to the loops ex- 
pansion when a magnetic field is applied to the sys- 
tem. Since quantum fluctuations favor collinear spins, 
one generically expects the magnetization to field curve 
to include plateaus at certain rational values of the 
magnetization, corresponding to collinear spin arrangc- 
ments^LiSiiS. Thus, at certain values, the field will 
choose collincar states satisfying 



Vi 



M , V simplexes a . 



(6.3) 



for some non-zero M. These are the new classical ground 
states in this theory, and now the question we ask is which 
of these are selected by harmonic quantum fluctuations. 
Repeating the derivation of the loop expansion l|4.4|l . wc 
now find that the diagonal elements of H change, as now 
Eq. lt2~TIjl is modified to H = ±W T W - r]M. Writing 
down the equations of motion in terms of the diamond 
lattice, as before, we find that Eq. I|4.1|) is still valid, but 
that the diagonal elements of fi are all non-zero and equal 
to —Al. In order to remove these diagonal elements, we 
can define 



fx = (j, + Ml , 



(6.4) 



where fi is equal to the M — value of fi as defined in 
Eq. H4.2f) . and, as before, it only connects nearest neigh- 
bor tetrahedra. The loop expansion is now 

-SN S + \M\N Z S 
= SVA^CaTr A^-2M^,) + ^-4A)l 



n=0 



-SN S + \M\N Z S 

oc ^ y n n — k / 
n=0 fe=0 j=0 v ' 

x(M 2 - 4A) n ~ k - j {-2M) jr TTH 2 k+3 
-SN S + \M\N Z S. 



(6. 



Here N z is the number of zero modes of W T W, i.e. N s /2 
for the pyrochlore lattice. Now the calculation goes as in 
Sec. IIV Al noting that the trace is non-zero only for even 
j. When one re-sums the terms in Eq. (|6.5|l to construct 
an effective Hamiltonian of the form H4.6f) . one finds that 
unlike the M = case, where always sign(if2z) = (— f )' +1 , 
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FIG. 11: (Color online) (001) projection of a pyrochlore lattice 
M = 2 0-fiux ground state. Each tetrahedron has 3 up spins (filled 
circles) and one down spin (empty circle). This is the ground state 
with the smallest possible (4 site) magnetic unit cell. The dark 
lines show x-z planes (coming out of the page), that can each by 
flipped independently to obtain a valid 0-flux M = 2 state. Since 
there are O(L) such planes, the ground state entropy is at least of 
order L. 



the signs of the expansion terms are no longer easy to 
predict. 

Applying this calculation to the only non-trivial 
collinear case on the pyrochlore lattice, i.e., M = 2, and 
rearranging the terms in the form 1)6. 2J) . we find (see 
Tab. [l| that B < 0, and that the coefficient of the inter- 
action terms, J , is two orders of magnitude smaller than 
the in effective field B. Therefore the complementary lat- 
tice ground state is a uniform fj a = 1 state, corresponding 
to the family of states with positive hexagon products 
(the so-called zero-flux states). 

As in the zero magnetization case, studied in Sccs.ll"fl- 
IfVI the harmonic energy is invariant under any gauge 
transformation that flips the spins in a set of tetrahc- 
dra, while conserving the constraint (|6.3|) . However, it is 
clearly much harder, in this case, to construct a transfor- 
mation in this way, because a single tetrahedron flip vio- 
lates the constraint not just on the neighboring tctrahc- 
dra, but on the flipped tetrahedron itself as well. There- 
fore one would expect the ground state degeneracy for 
this model to be smaller than for the M = case. 

It is easy to show that still, the ground state entropy 
is at least of order L, by observing that in the simplest 
ground state, with a 4 site unit cell, on can construct 
gauge transformations by independently flipping some 
of O(L) parallel planes, each composed e.g. of paral- 
lel antiferromagnetic lines in the x-z plane (see Fig. Illjl . 
The keys to these transformations being valid are: (i) 
The spin-product around any hexagon (and therefore any 
loop) in the lattice is not affected by it. (ii) Since an- 
tiferromagnetic lines arc flipped, the spin sum on each 
tetrahedron remains the same. 



B. Other lattices 

The calculation of the loop expansion is actually quite 
general and can be applied on any lattice composed of 
corner-sharing simplexes, as long as the classical ground 
state is collinear, and the (zero field) Hamiltonian has 
the form (|2.14|l . We have implicitly assumed that the 
simplex lattice, formed by centers of simplexes, is bipar- 
tite, although we could easily modify the calculation to 
take care of a more general case. Given these assump- 
tions, the only lattice information relevant to our calcula- 
tions is the coordination z of the simplex lattice and the 
lengths of the various loops in the same lattice. In most 
cases, the coordination of the simplex lattice is equal to 
the number of sites in a single simplex, e.g. z = 4 for the 
pyrochlore, or 3 for the kagome. However the expansion 
works equally well for cases where not every lattice site is 
shared by two simplexes, as in the capped kagome model 
below, in which case z is smaller than the number of sites 
in a simplex. 

In general, the ratio \K2i+2 / K2i\ is expected to increase 
and becomes closer to 1 for smaller values of z, meaning 
that we must take longer and longer loops into considera- 
tion, to determine the ground state. We find numerically 
that the convergence rate of the expansion l|4.4[) decreases 
as z decreases, as does the value of A at which conver- 
gence is obtained. The accuracy of our Bethe lattice 
approximation is expected to be better for larger values 
of z, since there are relatively fewer uncounted paths. 

1. Checkerboard lattice 

The checkerboard lattice is also often called the pla- 
nar pyrochlore, and is a two-dimensional projection of 
the pyrochlore lattic e) 34 ! 35 . It is composed of "tetrahc- 
dra" (crossed squares) whose centers form a square lat- 
tice (see, e.g. the (001) projections of the pyrochlore in 
Figs. |^|a, |3l4l9|) . Aside from the dimensionality, a major 
difference between the checkerboard and the pyrochlore 
lattice is that in the checkerboard, not all of the tetrahe- 
dron bonds are of the same length, and thus one would 
generally expect coupling J 1 ^ J along the square diag- 
onals. Nevertheless, to pursue the frustrated analog of 
the pyrochlore, we shall consider here the case J' = J. 

As in the diamond lattice, the coordination of the 
(square) simplex lattice is z = 4, and therefore the calcu- 
lation of approximate loop expansion coefficients, is iden- 
tical to the pyrochlore calculation of Sec. lIV Al However, 
we should note that since the shortest loops are now of 
length 4, the error in our Bethe lattice approximation is 
greater than in the pyrochlore case, since the "repeated 
loops" that we ignore carry more significant weight. Nev- 
ertheless, the ignored terms are still expected to be two 
orders of magnitude smaller than the Bethe lattice terms. 

In the checkerboard case, the complementary lattice 
is a square lattice composed of the centers of the empty 
square plaquettes. Now the effective field in Eq. i|(j.2fl is 
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the 4- loop coefficient B = K4/ S <0, which prefers ?7a>0, 
or in other words, 0-flux order. While the nearest neigh- 
bor complementary lattice coupling is antiferromagnctic 
J = Kq/S > and competes with B, it is not strong 
enough (see Tab.[I} to frustrate the uniform 0-flux order. 

The ground state entropy of the checkerboard lattice, 
has been shown in Ref. [la to be of order L, by construc- 
tion of all of the possible even and odd gauge transfor- 
mations on a particular reference 0-flux state. A simple 
explanation for the degeneracy is that, given a line of 
spins, say in the x direction, there are, at most, 2 choices 
in the construction of an adjacent parallel line. This is 
because there is one constraint on each tetrahedron (2 
down-spins) and one on each plaquette (an even number 
of down-spins). Thus, starting from an arbitrary choice 
of one horizontal line (2 L choices) , there are < 2 L ways of 
constructing the rest of the lattice. The resultant entropy 
is 0(L). 

Applying a magnetic field that induces the M = 2 
plateau, we find (Tab. [Q), that the complementary lat- 
tice effective Hamiltonian has B > 0, J < 0, favoring the 
7r-flux uniform state. The ground state entropy in this 
case is also of order L. To show this, we note that there is 
only one down-spin in each tetrahedron when M = 2, and 
that the spin product around each plaquette is —1. This 
means that there must be exactly one down-spin around 
each plaquette. 

We could use an argument similar to the one used in 
the M = model, to find the entropy in this model. A 
more elegant argument uses a 1-to-l correspondence be- 
tween the ground states and complete tiling the checker- 
board lattice with squares of size 2a x 2a, where a is the 
size of each plaquette (and "tetrahedron"). Here each 
square is centered on a down-spin and covers the two 
plaquettes and two tetrahedra to which it belongs. The 
entropy of such tilings is clearly of order L. 



2. "Capped kagome" 

The kagome lattice Heisenberg model has been one of 
the most studied highly frustrated models. The lattice is 
two dimensional, and is composed of corner sharing tri- 
angles, such that the centers of the triangles form a hon- 
eycomb simplex lattice. This model too is closely related 
to the pyrochlore, as a (111) projection of the pyrochlore 
lattice contains kagome planes sandwiched between trian- 
gular planes, such that the triangular lattice sites "cap" 
the kagome triangles to form tetrahedra. 

We cannot apply our collincar loop expansion to the 
kagome Heisenberg model, with no applied field, be- 
cause there are no collinear states that can satisfy the 
zero triangle-sum rule. One way to consider a "collincar 
kagome" is to look at a capped kagome lattice which con- 
sists of a kagome, flanked by two triangular lattices, so 
that each triangle turns into a tetrahedron. This model 
was studied by Tchcrnyshyov et. al^, who referred to 
it as a "[111] slice of pyrochlore". Those authors found 



that the ground state is one in which one out of every 
4 hexagons has a positive spin product (the 2x2 state 
shown in Fig. 112b). Surprisingly, the Hamiltonian for 
this model can be written in the form (|2.14|) . and there- 
fore we expect our loop expansion to work. Furthermore, 
as long as there is no applied field, B = Kq/S is positive, 
so applying our intuition based on the previously dis- 
cussed models, we would naively think that the ground 
state should be a 7r-flux state. 

However, this is the one model that we have studied, 
where the complementary lattice effective Hamiltonian 
has interaction terms strong enough to resist the B field 
term. Here, the complementary lattice nearest neighbor 
term corresponds to loops of length 10 and therefore J> 
(there are no loops of length 8 in the kagome), and 
J' = K 12 /S<0. 

In the 7r-flux state, all complementary lattice spins are 
— 1, and thus while the B term in Eq. (|6.2() is optimized, 
all of the nearest neighbor bonds (3A^, which N° is the 
number of complementary lattice sites) are violated, as 
well as all of the 3-site terms (2N^). On the other hand, 
the 2x2 state has only iV^/4 negative spins, but half of 
the complementary lattice bonds and 3/4 of the 3-spin 
terms are satisfied. Applying the coefficients we obtained 
(from Tab. [J J /B ~ 0.14), and including also the next 
order (K12) term JzjB « —0.05) to the expansion, we 
find the energy per complementary-lattice site 

— flux _ 77*2X2 

harm harm _ _q + 3 j _ 3^ _ 0.0026 > , 

(6.6) 

and we find that the 2x2 state is favored over the 7r-flux 
stated. We provide this calculation as an illustration 
for the difficulty in using of the effective Hamiltonian 
in a model with an unusually frustrated complementary 
lattice. In order to actually determine the ground state in 
this case, one must include further terms in the effective 
Hamiltonian. In fact, based only on the terms included 
in Eq. (|6.6|) . we would conclude, mistakenly, that the 
complementary lattice ground state is the \/3 x -\/3 state 
(not to be confused with the well-known coplanar kagome 
ground state with the same ordering vector), where 1/3 
of the spins are up, so that each triangular plaquette has 
2 down spins and one up spin (see Fig. 112b ). However, 
based on numerical diagonalization, we find that both 
the 2x2 state and the 7r-flux state have, in fact, lower 
energy than the VSx^/S state, in agreement with Ref. 

Due to the large number of degrees of freedom in this 
model, arising from the "free" spins capping each trian- 
gle, there is an extensive number of 2 x 2 ground states, 
as has been calculated in Ref. HBl 



3. Kagome with applied field 

Another way of obtaining a "collinear kagome" model, 
is to apply a field strong enough to induce collinear states 
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FIG. 12: (Color online) (a) The 2x2 ground state of the "capped 
kagome" lattice, (b) The V3 X V3 state in the same lattice. Up 
(down) spins are marked by filled (empty) circles, and ferromag- 
netic (antiferromagnetic) lattice bonds are shown in dark (light). 
The dashed lines represent the complementary lattice bonds, whose 
sites are marked + and — . 




FIG. 13: (Color online) The M = 1 kagome model ground 
states can be mapped 1-to-l to continuous lines on triangular 
lattice, such that a line can turn, at each site by ±60° or 
0°. The "kinks" in each lines must lie along directed defect 
lines (dashed lines). There are two possible cases: (a) parallel 
defects alternating in directions, (b) three lines meeting at a 
point defect (and no other kinks). 



with M = 1 on the (standard) kagome lattic o 33 i 36 . Ap- 
plying our expansion, we find B < and J < (see 
Tab.QJ, consistent with a uniform zero- flux ground state, 
as we have indeed confirmed by numerically calculating 
the zero-point energy. 

In the M = 1 kagome model, 0-flux states are obtained 
by spin arrangements in which there are exactly 2 down- 
spins around each hexagon. To find the ground state 
entropy, we map these to a dimer covering of the (trian- 
gular) complementary lattice, in which there are exactly 
2 dimers touching each complementary lattice site, and 
there is exactly 1 dimer in each plaquette. These can 
be viewed as continuous lines on the triangular lattice, 
that can turn at each node by ±60° or 0°. Since there 
is no way of closing such lines into loops, one finds that 
they generally run in parallel, possibly turning by 60° 
to form a "directed line defect" . Three such line defects 
can meet and terminate at a point, as long as they form 
120° angles and are all directed into this "point defect" . 
Thus, each ground state either has exactly one point de- 
fect and three line defects coming into it (and no others), 
as in Fig. 113b . or parallel line defects (alternating in di- 
rection), as in Fig. 113b . allowing for entropy of order L. 
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TABLE I: Effective complementary lattice Ising Hamilto- 
nian 16.21 coefficients for the various models that we consider. 
The coefficients were obtained by the expansion 16.511 . with 
n > 30 and extrapolated to n — > do. The complementary lattice 
ground state is given in the right column, and it corresponds, 
in all cases, to a family of states in the direct lattice. 



VII. CONCLUSIONS 

In this paper, we have presented a thorough semiclas- 
sical linear spin- wave theory of the pyrochlore lattice an- 
tiferromagnetic Heisenberg model. We found that for 
any collincar classical ground state, half of the spin-wave 
modes are generic zero modes, i.e. modes that have no 
dispersion, and that are the same for any state. Of the 
remaining N s /2 ordinary modes, some are optical in na- 
ture, while others vanish along the major axes in the 
Brillouin zone. These non-generic zero modes carry di- 
vergent fluctuations. We find that the divergent modes 
are closely related to gauge-like transformations that re- 
late exactly divergent (to this order in l/S) states. 

In Sec. II VI we studied the (harmonic) zero-point en- 
ergy and find that it can be expanded in loops, to ob- 
tain an effective Hamiltonian (|4.6|) written in terms of 
flux variables. We find the effective Hamiltonian coeffi- 
cients using re-summation of the Taylor expansion, via 
a Bcthe lattice approximation, in which we take advan- 
tage of the fact that all state-independent terms can be 
evaluated to excellent accuracy using closed paths on a 
coordination-4 Bethe lattice. We have calculated the har- 
monic zero-point energy using numerical diagonalization 
of the spin-wave spectrum to demonstrate that the ef- 
fective Hamiltonian provides us in good agreement with 
exact results. 

In Sec. [V] we showed that any non-collinear deviation 
away from a collinear state results in an energy increase, 
and thus justified the search for a ground state among 
collincar states only. While we did show that the collincar 
states are local energy minima in the space of all classical 
ground states, we did not show that there are no other, 
non-collinear minima. More work would be required to 
rigorously prove this point. 

It is instructive to compare our results to the well- 
known results on the kagome lattice. In both cases, there 
is an extensive classical degeneracy, associated with all 
states satisfying Eq. i|1.3fl . Whereas in the pyrochlore 
collincar states are favored by the zero-point energy, in 
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the kagome, which is composed of triangles, there are 
no collinear classical ground states. It turns out from 
the harmonic spin- wave theory of the kagome, that all 
coplanar states satisfying (|1.3|) are exactly degenerated, 
so that the entropy remains extensive. Here, on the other 
hand, we find that, as suggested in Ref.^J the pyrochlore 
harmonic order ground states are only a small subset of 
all of the collinear classical ground states, and that the 
C(L 3 ) classical entropy is reduced to an entropy of, at 
most, O(LlnL) (and at least O(L)). 

In order to search for a unique ground state one must 
go to higher order terms in expansion H1.4|^ i. as has been 
done previously for the kagome lattice^. In that case, 
anharmonic corrections to the zero-point energy resulted 
in a selection of a unique kagome ground state. We ex- 
pect that the spin-wave modes with divergent Gaussian 
fluctuations described in this paper will play a significant 
role in the anharmonic ground state selection. 

In Sec. I VII we showed that the loop expansion approach 
to the effective Hamiltonian can be easily applied to the 
case of a magnetization plateau induced by a large mag- 
netic field, as well as to other bisimplex lattices that 
support collinear spin configurations as classical ground 
states. We applied this to find the effective Hamiltonian 
coefficients and ground state of the M = 2 pyrochlore, 
the M = and M = 2 checkerboard model, the M = 1 
kagome model, and the capped kagome models, which is 
the only case that we have encountered where the ground 
state fluxes are non-uniform. Together with Prashant 
Sharma, we have recentljsi^ employed a similar effective 
Hamiltonian on the large- iV Heiscnbcrg model on the py- 
rochlore. We believe that this type of approach can and 
should be used more often in determining the ground 
state of highly frustrated models with many degrees of 
freedom. 
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APPENDIX A: BETHE LATTICE PATHS 

Consider a Bethe lattice of coordination z, and Nb 
sites. We shall assume that Nb is infinite so that the 
translational symmetry of the lattice is conserved. We 
would like to find the number of paths of length k that 
start and end at a particular site. Since all sites are 
equivalent, the total number of paths would just be Nb 
times this quantity. 

Define the following values: 

• fk- The number of paths of length 2k starting and 
ending at a particular site a. 



• gk- Same as fk, but counting only paths that do 
not return to the origin until the last step. 

• fk'- Same as fk, but the origin a only has z— 1 
nearest neighbors. 

• gk'- Same as gk, but the origin a only has z — 1 
nearest neighbors. 

The number fk can be calculated from {gi : i < k}. E.g. 



/3 = .93 + 2.91.92 + .9i • 



(Al) 



One way of calculating the coefficients of these expan- 
sions in by means of generating functions 



9kX 



F = 1 + G+G 2 



1 



k=l 

on the other hand 



1-G 



fk 



1 d k F 



k\ dx k 



(A2) 



(A3) 



x=0 



fk can be calculated from {gi : i < k} in an identical 
fashion. 

We calculate gk, gk and fk by induction: 

• ,9o = ,9o = fo = !• 

• Find g k = z/fe-i, gk = {z - l)h-i- 

• Calculate fk from {gi : i < k}. 

Note that we always obtain ~ (z — l) k , gk ~ (z — l) fc , 
and gk ~ z(z — l) fe , and the prefactor depends only on 
k. Once {gk} have been calculated we find {fk}, and 
calculate the Bethe lattice harmonic energy by Eq. 1)4. 4|) 
(where the Bethe lattice takes the place of the TVe-site 
simplex lattice) 



arm 

(Bethe) = SV^^C„^ 



(-4) 



i — k 



n=Q fc=0 



A k 



-SN, 



fkN B 
(A4) 



or similarly for a non-zero magnetization plateau using 
Eq. (|6.5|) . -Bharm (Bethe) is our approximation for the 
constant term Eq in the energy. 



APPENDIX B: CALCULATING LOOP 
COEFFICIENTS 

Here we provide the details of our calculation the co- 
efficient K21 in Eq. (|4.6|) . Consider the terms in the ex- 
pansion (|4.4|l that involve loops of length 21 and no other 
loops. In the Bethe lattice approximation, we assume 
that all of the paths involving a loop, can be viewed as 
a decorated loops, i.e., a loop with self- retracing paths 
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FIG. 14: Diagrammatic representation terms of summation of 
and /fc, which is used in the loop calculation of Apps.[E]and[X] The 

difference between the two terms is that whereas /j. is computed 
on a coordination z Bethe lattice, in ft^ the origin has only z — 1 
neighbors, and all other sites have z neighbors. 



We want to count the number of paths of length 2k in- 
volving a particular simple loop of length 21, and no other 
loops, such that each site along the loop may be an ori- 
gin of a self-retracing path (diagrammatically shown in 
Fig. I15fl . As explained in Sec. IIV A~3l in order to avoid 
double counting, we must consider trees whose origin has 
only z—1 nearest neighbors. Fortunately, we have already 
calculated such terms in App. lAl i.e. the terms /j. 

All we have to do is to find all of the ways of distribut- 
ing k — I steps that are not part of the loop, among 21 
sites, and take the product of / for each of those. In 




FIG. 15: Diagrammatic representation of the paths included in 
calculation of the hexagon coefficient. Each node along the loop is 
"dressed" by a Bethe lattice factor, fk, as shown in Fig. 1141 



(i.e. equivalent to Bethe lattice paths) emanating from 
each site. This means that within our approximation, all 
loops of length 21 are equivalent, and therefore, for any 
k, Tr fx, 2k in the expansion is equal to &21 multiplied by 
/[, the number of decorated paths of length 2k involving 
a particular loop of length 21. The effective Hamilto- 
nian (|4.6[) coefficient can therefore be written 



more concrete terms, for a given k, the number of possi- 
ble paths involving a particular loop of length 21 is 



n=0 k=0 




(Bl) 



In the case of a large magnetic field of Sec. lVl Al Eq. i|Bl|) 
would be replaced by the appropriate expression from 
Eq. jo3|) . 



A 




(k-l) 



fh ft2 

hi ii\ 



A, 



(B3) 



Cn 



i — k 



n 



K 2l (M) = s^E^EE, 

n=0 k=0 j=0 \ 



x (M 2 - AA) n ' k - j (-2M) j f k . (B2) 



We have multiplied the sum by 4/c because we can start 
anywhere along the path and go in any of two directions. 
Plugging the results, into Eq. I|B1(I (or l)B2|0 . we obtain 
the effective Hamiltonian coefficients. 
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